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摘要 炮 变 体 具有 明显 的 、 剧 烈 的 大 幅度 光 变 ， 中 长 时 标的 光 变 研究 对 于 揭示 泡 变 


体 的 光 变 特征 和 光 变 机 理 有 重要 作用 。 通 过 选取 平 谱 射电 类 星体 3C 454. 3 的 光学 B、V、 


R 和 红外 J、 波段 的 2008 年 6 月 ~2017 年 7 月 的 原始 光 变 数据 ， 采 用 功率 谱 方 法 ， 研 
究 3C 454.3 的 中 长 周期 光 变 特性 ， 研 究 得 出 光 变 主 周期 为 1. 25 年 ，4. 57 年 的 周期 为 
1.25 年 周期 的 厂 加 。3C 454.3 在 光学 、 红 外 以 及 射电 波段 的 光 变 有 一 定 的 关系 。 研 究 


显示 出 3C 454. 3 的 红外 光度 比 光学 波段 更 为 明亮 ， 红 外 光 变 比 光 学 波段 更 为 剧烈 。 


关键 词 ” 平 谱 射 电 类 星体 ，3C 454. 3， 长 周期 光 变 ， 功 率 谱 方 法 
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1 引言 


3C 454. 3( 红 移 z=0. 859) 是 一 个 低 光 度 、 低 红 移 ， 而 且 有 明显 剧烈 光 变 的 平 谱 射电 


类 星体 。3C 454. 3 是 多 个 小 组 都 很 关注 的 类 星体 ， 它 在 光学 波段 的 观测 已 经 有 了 很 长 的 


历史 ， 有 许多 小 组 对 3C 454. 3 进行 了 大 量 的 观测 和 理论 研究 。 有 研究 表明 3C 454. 3 的 


光学 波段 存在 越 红 越 亮 的 趋势 。3C 454. 3 在 伽 马 波 段 、 光 学 波段 和 红外 波段 


都 有 强烈 的 光 变 ”““， 其 在 红外 波段 的 剧烈 光 变 与 伽 马 射线 源 有 关 “。3C 454. 3 的 光 变 


曲线 存在 长 周期 光 变 ， 甚 至 有 的 耀 变 体 存在 十 多 年 的 长 周期 光 变 ””， 长 期 的 观测 数据 


对 研究 光 变 周期 及 其 辐射 模型 有 重要 意义 ””。 目 前 ， 功 率 谱 方法 常用 于 天 体 光 变 周期 


性 的 讨论 ， 张 卑 晶 中 等 人 曾经 用 功率 谱 方法 计算 BL Lac 天 体 S5 0716+714 的 B、V、R、 
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I 波段 数据 得 到 了 1. 1 年 的 长 时 标 光 变 。 Li 等 人 ”利用 Jurkevich 方法 和 功率 谱 方法 研 


究 了 3C 454.3 的 射电 22GHz 和 37GHz 波段 存在 (1. 57 土 0. 12) yr 及 (6.15+0.50) yr 


的 周期 。 光 变 研 究 对 于 揭示 炮 变 体 的 光 变 特征 和 光 变 机 理 有 重要 作用 。 
2 3C 454. 3 长 期 光 变 曲线 

通过 SMARTS 数据 库 获 取 了 2008 年 6 月 至 2017 年 7 月 3C454.3 的 光学 和 红外 波段 
接近 10 年 的 历史 光 变 数据 ， 其 中 包括 了 光学 B、V、R 和 红外 J. 波段 ， 分 别 有 891、 
855、878、861 和 752 个 观测 数据 ， 共 计 有 4237 个 观测 数据 。 如 果 每 天 均 可 以 观测 ， 则 
光学 BV、R 和 红外 JK 波 段 的 实际 观测 夜 占 理 论 观测 夜 的 履 盖 率 分 别 为 26. 81%、25. 72%、 


26. 41%、25. 90% 和 22. 62%， 平 均 为 25. 49%， 侠 有 如 此 多 波段 密集 观测 的 天 体 。 若 以 天 


为 采样 单位 ， 则 采样 频率 约 为 0. 25 次 /天 ,依据 Nyquist 采样 定理 ， 采 样 周 期 要 小 于 最 


小 信号 周期 的 一 半 ， 采 样 频率 应 该 大 于 最 大 信和 号 频率 的 两 倍 ，Nyquist 频率 约 为 0. 125 


次 /天 ， 可 以 得 知 ， 此 观测 数据 用 来 计算 该 天 体 存 在 大 于 8 天 的 光 变 周期 符合 Nyquist 


采样 定理 。 
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图 1 B、V、R、J、 下 波段 的 光 变 曲线 


Fig.1 Light curve in B, V, R, J and K bands 
图 1 为 10 年 间 3C 454. 3 AB. V. R J. KY BOAR HHA. TE 2008 年 6 月 至 2017 


年 7 月 期 间 ， 该 天 体 B 波段 的 星 等 平均 值 为 16. 108mag， 星 等 最 大 值 为 17. 337， 星 等 最 


小 值 为 13. 964, A 4473. 373mag; V 波段 的 星 等 平均 值 为 15. 537mag, 星 等 最 大 值 为 16. 756, 


星 等 最 小 值 为 13. 448，A as=3. 308mag; R 波段 的 星 等 平均 值 为 15. 122mag， 星 等 最 大 值 


为 16. 585， 星 等 最 小 值 为 12. 97， A as=3. 615mag; J 波段 的 星 等 平均 值 为 13. 419mag, 


T 星 等 最 大 值 为 15. 092， 星 等 最 小 值 为 11. 012， ^54-4.08mag; K 波段 的 星 等 平均 值 为 


X 15.017, 星 等 最 小 值 为 9. 188，A as=5. 829mag。B、V、R、J、 


Ex 


11. 594mag， 星 等 最 大 


K 波段 的 变化 量 至 少 有 3 NET. WK 1。 


表 1 B、V、R、J、 下 波段 星 等 变化 情况 


Tabel1 Variation of magnitude in B, V, R, J and K bands 


X Pt B V R J K 

c 最 大 值 (mag) 17. 337 16. 756 16. 585 15. 092 15. 017 
最 小 值 (mag) 13. 964 13. 448 12. 970 11. 012 9. 188 

平均 值 (mag) 16. 108 15. 537 15. 122 13. 419 11. 594 

变化 量 (mag) 3,919 3. 308 3. 615 4. 080 5. 829 
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图 2 B. V. R J. 人 波段 星 等 最 大 值 、 最 小 值 、 平 均值 和 差 值 


Fig. 2 Maximum, minimum, average and difference of magnitude in B, V, R, J and K bands 


在 光 变 曲线 中 ， 观 测 数据 平均 误差 约 为 0.02mag， 光 变 曲 线 测量 值 区 间 平 均 为 ^ es 
=4. 04lmag， 信 噪 比 为 S=0. 5%， 信 噪 比 低 ， 我 们 在 光 变 曲线 和 周期 分 析 中 未 考虑 观测 误 
差 所 带 来 的 影响 。 


从 光 变 曲线 可 以 看 出 3C 454. 3 天 体 的 光学 和 红外 波段 光度 在 十 年 间 变化 十 分 剧烈 。 


我 们 用 3C 454.3 的 B、V、R、J、 玉 波段 的 中 心 波长 来 研究 其 于 对 应 星 等 最 大 值 、 星 等 


最 小 值 、 平 均值 和 差 值 的 关系 。 根 据 图 2 来 看 ， 显 示 出 在 B、V、R、J、 上 波段 中 的 波长 


越 长 ， 对 应 的 星 等 最 大 值 、 星 等 最 小 值 和 平均 值 有 越 低 的 趋势 ， 然 而 对 应 的 差 值 ( 星 等 
变化 量 ) 有 越 大 的 趋势 。 图 1 中 还 显示 出 3C 454.3 WB V R J KEREC HR 
升降 走势 近乎 一 致 。 


3 3C 454. 3 的 光 变 性 质 


文中 采用 功率 谱 方法 讨论 3C 454. 3 的 周期 。 


3. 1 功率 谱 分 析 方法 


在 功率 谱 方 法 分 析 中 ， 若 存在 一 个 平稳 随机 信号 x( 的 自 协 方差 函数 为 Yx:(5)， 并 


满足 
太 17yx(z)ldzr < oo (19 


定义 x(b 的 功率 谱 密度 函数 为 yY(ZT) 的 傅 里 时 变换 ， 即 


Sole J- rdre dr (2) 
SS fe LY fu 
Y(T) == J", Sw) ei?'do. (3) 
由 自 协 方差 函数 的 定义 有 
Tx(T) = lima zy J ple) = helll — 1) = prl dr. (4) 


由 于 Sx(w)dw 表 示 在 @ 与 @ 十 do 之 间 信 和 号 所 提供 的 功率 。 故 Sx(w) 代 表 功 率 相对 与 


频率 的 分 布 ， 则 称 为 功率 谱 函 数 密度 。 
对 于 平稳 时 间 序 列 {Xn}, 如 它 的 协 方差 序列 (Kk), k=0, 1, 2,… 满 足 


È=- |Y XK)| < oo (5) 


Jus SCOGMT] EL SET e SE AJ 
miuw)e*x Yk)e en, (6) 


上 上 式 中 At 为 数据 采样 间隔; 上 自 协 方差 序列 XYx(k) 定 义 为 


Y (k) = lim, ,, ee (7) 


功率 谱 密度 函数 也 可 以 直接 从 数字 序列 的 离散 傅 里 时 XCD 利 用 


(w) = GX)? (8) 
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HA | Fi] | 


获得 。 


因数 字 序 列 的 离散 傅 里 叶 变 换 具有 周期 性 , c (oo) TAA AWE, PARWA 


数 季 序 列 的 周期 图 估计 。 


可 证 ""， 对 于 一 个 随机 信号 


Y, genes + an (9) 


an 为 平稳 正 态 序列 , 其 功率 密度 Sa(w) 在 oo 处 连续 , PTVAY, B3 58] SHAM, (w) de [n m] 


上 oo 附近 的 某 一 频率 处 达到 最 大 值 ， 并 且 当 WW 很 大 时 有 


00 , W = Op» 
jim In(W) = E (10) 
0,0209. 


(10) 式 中 , wo 是 信号 ;的 固有 频率 。 当 NRK, n (WEO MEER- SSE e 


yc 在 其 他 的 频率 处 功率 谱 的 密度 值 较 小 。 根 据 这 种 特征 ， 就 可 以 寻找 信号 中 的 周期 成 分 。 


当 信 号 含有 不 止 一 个 周期 分 量 时 ，h(w) 中 会 出 现 多 个 极 大 值 。 为 了 判别 真 伪 周期 ， 则 


a 


N 需要 根据 周期 图 极 大 对 应 的 频率 值 , 利用 最 小 二 乘法 求 出 相对 应 各 频率 的 周期 分 量 的 振 


幅 来 比较 大 小 。 


SMARTS 数据 中 的 采样 数据 并 非 均匀 采样 时 间 序 列 ， 用 功率 谱 方法 计算 3C 454.3 的 


光 变 周期 时 对 光 变 数据 的 空白 处 进行 了 相 邻 数据 平均 值 插值 处 理 。 


以 上 为 功率 谱 方法 测量 周期 的 基本 原理 ， 其 对 于 耀 变 体 长 周期 光 变 研究 可 靠 性 的 理 


论 见 论文 [9]。 


3.2 结果 


通过 功率 谱 方法 计算 ， 得 到 了 3C 454.3 fj B. V. R J. K 波段 功率 谱 图 ， 图 中 显 


示 了 2 个 主要 的 频率 峰值 。 
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从 B 波段 功率 谱 图 来 看 ， 存 在 周期 有 454.5 和 1666.7 天 ， 分 别 对 应 1.25 年 和 4. 57 
年 。 峰 值 频 率 为 0.00060 (1/day) 的 半 高 宽 为 0.00035; 峰值 频率 为 0. 00220 (1/day) 


的 半 高 宽 为 0.00030。 如 图 3 所 示 。 


iy 4 
a | l 1 | 1 | 


i l 
0.005 0.006 0.007 0.008 0.009 0.01 
Frequency(l/day) 


图 3 B 波段 功率 谱 图 


Fig.3 The power spectrum figure in B band 


从 V 波段 功率 谱 图 来 看 , 存在 周期 有 454. 5 和 1666. 7 R, 分 别 对 应 1. 25 年 和 4. 57 
年 。 峰 值 频 率 为 0.00060 (1/day) 的 半 高 宽 为 0.00038; 峰值 频率 为 0. 00220 (1/day) 


的 半 高 宽 为 0.00032。 如 图 4 所 示 。 


1 | 1 | l | l | 
0 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 
Frequency(l/day) 


图 4V 波段 功率 谱 图 
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Fig.4 The power spectrum figure in V band 


从 RR 波段 功率 谱 图 来 看 , 存在 周期 有 454. 5 和 1666.7 天, 分 别 对 应 1. 25 年 和 4. 57 
年 。 峰 值 频 率 为 0.00060 (1/day) 的 半 高 宽 为 0.00035; 峰值 频率 为 0. 00220 (1/day) 


的 半 高 宽 为 0. 00030。 如 图 5 所 示 。 


| 1 | 
0 0.001 0.002 0.003 0.004 0.006 0.007 0.008 0.009 0.01 


0.005 
Frequency(/day) 


图 5 R 波 段 功率 谱 图 


Fig.5 The power spectrum figure in R band 


从 了 丁 波段 功率 谱 图 来 看 ,存在 周期 有 454.5 和 1666. 7 R, 分 别 对 应 1. 25 年 和 4. 57 
年 。 峰 值 频 率 为 0.00060 (1/day) 的 半 高 宽 为 0.00036; 峰值 频率 为 0. 00220 (1/day) 


的 半 高 宽 为 0.00030。 如 图 6 所 示 。 
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图 6 J 波段 功率 谱 图 


Fig.6 The power spectrum figure in J band 


从 下 波段 功率 谱 图 来 看 ,存在 周期 有 434.8 天 和 1666.7 天, 分别 对 应 1.19 年 和 4.57 


年 。 峰 值 频 率 为 0.00060 (1/day) 的 半 高 宽 为 0.00038; 峰值 频率 为 0.00230 (1/day) 


的 半 高 宽 为 0.00033。 如 图 7 所 示 。 


1 | 1 | 1 | 


L l 
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图 7 波段 功率 谱 图 


Fig.7 The power spectrum figure in K band 


各 个 波段 的 功率 谱 图 峰值 频率 及 其 半 峰 全 宽 依据 表 2 所 示 。 通过 计算 得 到 3C 454. 3 


在 光学 和 红外 波段 都 存在 2 个 光 变 周期 ， 分 别 为 1.25 年 和 4.57 年 。 


表 2 B、V、R、J、 下 波段 的 功率 谱 图 峰值 频率 及 其 半 高 宽 


Tabel2 Peak frequency and FWHM of Power Spectrum in B, V, R, J and K bands 


波段 B V R J K 
峰值 频率 1 0. 00060 0. 00060 0. 00060 0. 00060 0. 00060 
峰值 频率 1 半 高 宽 0. 00035 0. 00038 0. 00035 0. 00036 0. 00038 
峰值 频率 2 0. 00220 0. 00220 0. 00220 0. 00220 0. 00230 
峰值 频率 2 半 高 宽 0. 00030 0. 00032 0. 00030 0. 00030 0. 00033 


4 讨论 


从 图 1 显示 出 ， 在 3C 454.3 AB. V. R J. K，5 个 波段 中 ， 历 史 光 变 曲 线 中 变化 


近乎 一 致 ， 图 3 一 7 显示 出 功率 谱 图 大 致 相同 ， 这 有 可 能 产生 于 3C 454. 3 在 光学 和 红外 


波段 的 光度 的 变化 根本 上 有 着 一 定 的 关系 。 


根据 图 2 来 看 ， 表 现 出 在 B、V、R、J、 下 波段 中 的 波长 越 长 ， 对 应 的 星 等 最 大 值 、 


x 


星 等 最 小 值 和 平均 值 有 越 低 的 趋势 ， 然 而 对 应 的 差 值 有 变 大 的 走势 。 星 等 越 低 ， 其 光度 


也 就 越 亮 ， 这 暗示 出 在 大 部 分 时 间 上 ， 红 外 波段 比 光 学 波段 更 为 明亮 。3C 454. 3 光度 变 


化 时 ， 红 外 波段 比 光学 波段 较为 


剧烈 。 Fan 错误 ! 未 找到 引用 源 。 错 误 ! 未 找到 引用 源 。 等 人 研究 了 3C 454,3 


在 2006— 2011 年 光学 波段 明显 存在 越 红 越 亮 的 趋势 , 这 个 趋势 到 红外 本 波段 也 适用 。 


依据 表 2 所 示 ， 各 波段 都 显示 出 2 个 显著 周期 ,频率 为 0. 00060 (1/day) 的 周期 在 


5 个 波段 是 完全 一 致 的 。 而 和 波段 434. 8 的 周期 不 同 于 其 余 4 个 波段 ， 通 过 计算 并 查 表 


可 得 , 各 波段 该 周期 均值 为 450.6 天 ,周期 置信 和 度 为 0.95 的 置信 区 间 为 [442.9，458.3] 天 ， 
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即 450.6 土 7.72 天 ， 区 间 可 表述 为 [1.21，1.26] 年 。 说 明了 1.25 年 的 周期 具有 较 高 的 置 


言 水 平 ， 可 以 作为 3C 454. 3 的 主要 周期 。 


一 个 


半 高 全 宽 较 小 的 峰值 对 应 的 误差 相对 来 说 较 小 ， 因 此 峰值 频率 为 0. 0022. C1/day) 


的 周期 比 峰值 频率 为 0.0006 (1/day) 的 周期 更 为 显著 ; 前 者 对 应 的 周期 为 1. 25 年 ， 后 
者 对 应 的 周期 为 4. 57 年 。1. 25 年 的 周期 与 Gorshkov" "等 人 在 射电 波段 21. 7GHz 计算 
得 到 的 1. 2 年 周期 光 变 大 致 相符 合 ,因此 1. 25 年 的 周期 可 以 辨认 为 3C 454. 3 的 显 周期 ， 


其 光 变 曲线 主要 由 1. 25 年 的 周期 导致 。 在 光学 波段 、 红 外 波段 和 射电 波段 的 周期 相同 ， 
可 以 认为 3C 454. 3 的 光度 变化 及 其 周期 特征 现象 在 光学 、 红 外 和 射电 波段 具 一 定 的 关 


ZINO 


本 文 使 用 功率 谱 方 法 得 到 的 周期 有 1. 25 年 和 4.57 年 ， 其 中 4.57 年 的 周期 因为 受 


观测 时 间 和 观测 数据 的 影响 ， 其 半 高 宽 略 大 ， 误 差 相 对 较 大 ， 与 其 真实 周期 可 能 会 有 偶 
差 ， 视 为 1. 25 年 周期 的 车 加 。 
5 结论 


对 3C 454. 3 在 B、V、R、J、 波 段 的 光 变 比较 分 析 ， 显 示 出 3C 454. 3 的 红外 波段 


的 光度 比 光学 波段 更 为 明亮 ， 红 外 波段 的 光度 变化 比 光 学 波段 更 为 剧烈 。 


我 们 用 功率 谱 方 法 对 FSRA 3C 454. 3 中 B、V、R、J、 波段 的 光 变 曲线 进行 研究 ， 


结果 表明 该 天 体 在 光学 、 红 外 波段 存在 光 变 主 周期 为 1.25 年 ，4. 57 年 的 周期 为 1.25 


年 周期 的 登 加 。 由 于 相对 以 前 的 研究 中 ， 信 和 号 采样 是 中 长 周期 光 变 的 研究 障碍 而 不 可 和 避 


免 ， 本 文 在 利用 最 新 (2008—2017 年 ) SMARTS 计划 的 观测 基础 之 上 ， 首 次 利用 功率 谱 


方法 计算 3C 454.3 的 光学 红外 B、V、R、J、 下 波段 周期 ， 获 得 了 相对 准确 的 结果 。 在 


光学 波段 、 红 外 波段 和 射电 波段 的 周期 相同 ， 可 以 认为 3C 454. 3 的 光度 变化 及 其 周期 
特征 现象 在 光学 、 红 外 和 射电 波段 具有 一 定 的 关系 。 
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The medium and long period light variation characteristics of 
FSRQ 3C 454. 3 
Wu YueCheng, Zhang HaoJing, Yu Lian, Xu Xiaolin 


(College of Physics and Electronics, Yunnan Normal University, Kunming 650500, China) 


Abstract Blazars have obvious and intense large-scale light variations. 
The study of light variations on medium and long time scales plays an 
important role in revealing the characteristics and mechanism of light 
variations of Blazars. By selecting the original optical B, V, R and 
infrared J, K bands of the flat-spectrum radio quasar 3C 454.3 from June 
2008 to July 2017, the medium and long period characteristics of 3C 454. 3 
are studied by power spectrum method. It is concluded that the main period 
of light variation is 1.25 years, and the period of 4.57 years is the 
superposition of 1.25 years. 3C 454.3 is related to the optical, infrared 
and radio variations. The results show that the infrared luminosity of 
3C 454.3 is brighter than that of optical band, and the infrared light 
variation is more intense than that of optical band. 


Key words FSRQs, 3C 454. 3, long-term variability, method of power spectrum 


